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I. INTRODUCTION 



The understanding of the kinematics and dynamics of the impulsively 
started flow about bluff bodies in general and a circular cylinder in 
particular is of both practical and fundamental importance. From a 
theoretical point of view, analytical and experimental results enable 
one to understand the evolution of the separation of flow and the sym- 
metric vortex development. 

From a practical point of view, the understanding of the character- 
istics of impulsive flow enables one to calculate the normal and side 
forces acting on slender bodies of revolution at large angles of attack 
through the use of the so-called "impulsive flow analogy". 

The present investigation addressed itself to the experimental 
determination of the drag (normal) force and lift (side) force acting 
on smooth cylinders immersed in an impulsively started flow through the 
use of a vertical water tunnel. 
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II. PREVIOUS STUDIES AND THE 
IMPULSIVE FLOW ANALOGY 



When a slender body is inclined to a flow, it experiences a compo- 
nent of force in the plane in which the inclination occurs and a com- 
ponent perpendicular to that plane. These forces will be referred to 
hereafter as drag and lift forces. Earlier works were mainly concerned 
with the drag forces at moderate angles of attack. The flow pattern at 
small to moderate angles of attack consists of a pair of axially 
oriented symmetric vortices on the lee side of the lifting body. This 
physical flow field was approximated by the so-called NACA vortex model 
(See Fig. I). Munk CO was one of the first to study the problem and 
to introduce the cross-flow approach, which has been the basis of all 
drag-force prediction methods. Munk's work was primarily concerned 
with an ideal fluid and only relevant to small angles of inclination. 
Allen and Perkins £,23 considered slender bodies at large angles of in- 
clination, where the viscous effects and separation are important, and 
modelled the cross-flow development by analogy with the flow over an 
impulsively started cylinder. In other words, they have hypothesized 
that the evolution of the spiralling vortex sheets with di stance a long 
a slender body is similar to the evolution of symmetric vortices w i th 
ti me in an impulsively-started flow about a circular cylinder. This 
hypothesis turned out to be a very fruitful one and has been developed 
by later workers; for example Sarpkaya [33, Kelly [4], Hill [53, and 
Bryson [63 devised various approximate discrete-vortex models to study 
the evolution of vortices behind an impulsively started cylinder. 
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« Figure I. Symmetric vortex separation on a slender body of 
revolution at high angle of attack, NACA model. 
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These models were able to predict in an approximate way only the early 
stages of motion. These investigators were not concerned with the 
development of asymmetry in the wake and subsequent vortex shedding. 

Vortices peel off from missiles flying at high angles of attack 
and give rise to significant side or lift forces. In spite of the 
realization of the need for theoretical analysis and experimental data 
either for side forces acting on missiles or for lift forces acting on 
cylinders in impulsively started flow, there has been practically no 
data or systematic study. 

The development of flow about a cylinder started impulsively from ' 
rest was first studied by Blasius C^II. He has shown that the flow 
about the cylinder does not begin to separate until a distance of 
s = 0.351 c is covered by the flow. Blasius' work has later been ex- 
tended by Goldstein and Rosenhead CS], Gortler [?!]> Schuh ClOH, 

Watson Cll3» and Wundt Cl2D« They have demonstrated that the separa- 
tion depends on the shape of the body, the Reynolds number, and the 
particular time dependence of flow from the start of motion. Follow- 
ing the inception of separation at the rear stagnation point, the 
separation points move rapidly around the body until at large values 
of time they reach a nearly steady state position, often quite close 
to that observed experimentally in steady flow about the same body. 

In spite of its relevance to the impulsive-flow analogy, relative- 
ly little experimental work has been carried out on impulsively 
started flow about cylinders. Schwabe Cl3D dragged a cylinder of 9 cm 
diameter in a shallow channel at a speed of 0.8 cm/sec through the use 
of weights and pulleys and attained Reynolds numbers of about 600. He 
developed a procedure which enabled him to calculate pressure 
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distribution and hence the drag force from the pictures of the flow 
pattern. Schwabe’s drag coefficient is about twice as large as the 
steady state value and is still increasing after the cylinder has 
displaced a distance of about 9 body radii. Approximate numerical 
analysis carried out by Bryson shows a maximum of about 1.5 for 
the drag coefficient at the time when the cylinder moved about 4 body 
radii. Subsequently, the drag coefficient drops to an unrealistically 
low value as noted by Bryson. 

The only other investigation in recent times was undertaken by 
Sarpkaya C3D through the use of a vertical water tunnel. He found 
that the drag coefficient reaches a value of about 1.5 at a relative 
fluid displacement of about 8 body radii. In spite of the difficul- 
ties encountered by him in generating an impulsively-started flow, 
his results have clearly shown that the drag coefficient does reach 
a maximum of about 1.5 and then rapidly decreases to a steady state 
value of about 1.2. The method developed by him has been further 
improved and utilized in the present investigation not only for the 
purpose of more accurately determining the exact position of the 
maximum drag coefficient but also the evolution of the lift coeffi- 
cient. Consequently, the work described herein presents drag force 
data which may be compared with those obtained previously and also 
lift force data for which there is no previous work. 
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III. EXPERIMENTAL ARRANGEMENT 



A. TEST FACILITY 

The experiments described herein were conducted in a vertical 
water tunnel designed by Professor T. Sarpkaya (See Fig. 2). A 
quick release valve mechanism which opened or closed with the help 
of a three-way valve was mounted at the bottom of the tunnel. When 
the quick release valve was opened, the fluid (water) moved down- 
ward through the tunnel into an underground reservoir. The partic- 
ular valve design is shown in Fig. 3. It allows an initial drop 
within a very short time period and then a slow further opening. 

The rate of initial drop as well as that of the subsequent valve 
motion is controlled by a piston in the lower part of the valve and 
the viscosity of oil in the liquid chamber. The amount of the 
initial drop, the area of the opening between the liquid chamber and 
the upper air chamber as well as the viscosity of the oil are vari- 
able and easily adjusted to obtain constant velocities at desired 
rates. Suffice it to note that the system was capable of producing 
an impulsively started steady flow with very little initial disturb- 
ances at desired velocities. In fact the acceleration of the flow 
was confined to a time period of about 0.05 sec for steady flow 
velocities of about 40 cm/sec. 

B. ACCELERATION AND VELXITY MEASUREMENTS 

The acceleration of the fluid was monitored by means of a differ- 
ential-pressure transducer connected to two pressure taps along a 
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Figure 2 
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Figure 3. Valve assembly. 
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vertical line and with a spacing The acceleration was calculated 
from the fact that the measured differential pressure Ap and the 
acceleration are related by 






dt 



Ap 



(I ) 



in which p is the density of water and represents the instantaneous 
acceleration. Repeated measurements of acceleration have shown that 
the velocity reaches a constant value within 0.05 sec and the accelera- 
tion remains zero thereafter. 

The velocity was measured through the use of a variable resistance 
probe. The change in resistance of an about 150 cm long platinum 
wire, placed vertically in the tunnel away from the tunnel walls, was 
recorded by an amplifier recorder assembly together with the drag or 
lift force. The slope of the elevation versus time curve yielded the 
velocity of the fluid motion through the use of the calibration data. 

C. TEST BODIES AND FORCE MEASUREMENTS 

Two smooth aluminum circular cylinders, 60.6 cm long with diameters 
62.23 mm and 75.82 mm, respectively, were used in the experiments. 

The cylinders were turned on a lathe from aluminum pipes. The length 
of each cylinder was such that it allowed a gap of about 2 mm between 
the tunnel wall and each end of the cylinder. The cylinder was pre- 
vented from moving towards one or the other wall by means of small 
0-rings attached to the round cantilever end of the force transducers. 

A double-ball precision bearing was inserted at each end of the 
cylinder in aluminum housings which sealed the cylinder air tight. 

The outer face of the bearing was flush with the cylinder. 
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Two Identical force transducers, one at each end of the cylinder, 
were used to measure the instantaneous drag and lift forces. A special 
housing was built for each gage so that the gage could be mounted on 
the tunnel window and rotated to measure either the drag or 'the lift 
force alone. The bellows protecting the strain gages were filled with 
Dow Corning 3140-RTV coating for water proofing. Then the ends of the 
bellows were sealed air tight. 

After mounting one of the cylinders, the exact angular position of 
each gage within its housing was determined by hanging a load on the 
cylinder and rotating the gage until the lift force component became 
zero. The final positions of the gages were marked with a pin. Finally, 
four bolts were placed on the gage housing to hold the gage rigidly 
in position. The removal of these bolts and the set pin allowed the 
rotation of the gage exactly 90 degrees to measure the drag force. 

Thus, both gages were capable of measuring either one or the other com- 
ponent of total resistance. Usually, one gage was set to measure the 
drag force and the other, the lift force. 

The calibration of the gages was accomplished by hanging loads at 
the mid length of the cylinder after setting both gages to sense only 
the drag component. The calibration was also checked by filling the 
tunnel with water and recording the signal generated by the buoyant 
force. This procedure provided an easy check on calibration and was 
used prior to each experiment. 
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IV. DISCUSSION OF RESULTS 



The analog data for both the lift and drag forces were evaluated 
through the use of the calibration constants and plotted in terms of 
the normalized distance Ut/c where U represents the velocity of the 
steady flow; t, the time from the start of the motion; and c, the 
radius of the cylinder. The drag and lift coefficients were defined 
as fol lows: 



and 



Drag Force 

T 1 

^ p L D U 



Lift Force 

1 ? 

^ p L D U 



( 2 ) 



(3) 



in which p represents the density of water; L, the length of the 
cylinder; and D, the diameter of the cylinder. 

A. DRAG COEFFICIENT 

Figures 4 and 8 show the drag coefficient as a function of Ut/c 
for various values of the Reynolds number defined by Re = UD/v. 

These figures show, regardless of the Reynolds number, that the drag 
coefficient always rises rapidly to a value of about I . 55 at Ut/c 

values from 4 to about 5 and then drops rapidly. The rate of drop as 

: \ 

well as the time at which the drag coefficient reaches a steady state 

/ ‘ 

value appears to depend on the Reynolds number. At relatively lower 
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Figure 4. C„ vs Ut/c for Reynolds numbers 15,350, 16,200, and 16,650. 



[ll 




2.001 



































i 

1 
























+ L 














r 
























t ^ 

+ t 




































4. 




































1 

+ , 




































+ ’ 
•k ' 














c 

e 

pr 


0 K) 
























r 












a s) 

- CO 






















R 


1 










CNJCNJ 
II II 






















> 

r > 












■ LU LU 

cc cn 
























r 














V *r 

1 


1 






















1 .X 


















r '“■' 

1 

1 






















+x 






L 
































t 








\ 
















— 
















1 




















1 

! j 


I 
















X 




















i 

1 




















y^\ 




















































































































































1 








y-rf^ 








































\! 
























— 








M 

i 


k 


“1" 


X 


X 


































+ 


+ 


■H 


' 

+■ ‘ 


< 

+ 


X 

} 




1 


X- 


, ___v .#1 








i 

— ( 


























+ 1 




^ -1 


A 




J 




*1 


i 
























1 


1 


1 

1 




1 

1 





§ 8 



cs ca 



fSI 



G3 

<\» 

(N 



C3 

• 

Q 

CM 



cs 

00 



ca 

CO 



CO 



ca 

rsl 



ca 

ca 



cs 

00 



C 3 

CO 



ca 



ts 

CM 



ca 



CJ 



(S I 



lN3I3IdJ303 OUyO 



1 



21 



Figure 5. C„ vs Ut/c for Reynolds numbers 20,300 and 20,700. 
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Figure 6. C vs Ut/c for Reynolds numbers 32,000, 32,600, and 33,200. 
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Figure 7. vs Ut/c for Reynolds numbers 32,900 and 32,800 (two runs). 
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Reynolds numbers (Re = 17,000), the drag coefficient drops to its 
steady state value at about Ut/c = 7 and thereafter fluctuates about 
C|^ = l.l as a result of the lift-induced drag force oscillations. At 
relatively higher Reynolds numbers (Re = 33,000), the maximum peak 
in the drag coefficient occurs at Ut/c = 5 ; and its value is slightly 
lower than that for the lower Reynolds numbers. One of the most sig- 
nificant differences between the evoiution of the drag coefficient 
for the smaller and larger Reynolds numbers Is in the manner whereby 
the drag coefficient reaches its steady state value. As noted earlier, 
at lower Reynolds numbers rapidly attains its steady state value 
at about Ut/c = 7. At higher Reynolds numbers, however, first de- 
creases from its peak value to a new plateau between 1.25 and 1.3 for 
Ut/c values from 7 to 13 and then decreases once more to its steady 
state value. 

The foregoing discussion covered only the general character! sties 
of the evolution of the drag coefficient with time without making an 
attempt to explain the reasons giving rise to them. The present ex- 
periments as well as those previously conducted by Sarpkaya C33 show 
that during the initial and symmetric development of vortices (as 
evidenced by the absence of lift force as will be discussed later) 
the vorticity rapidly accumulates in the wake. The symmetric vortices 
grow to sizes considerably larger than those which would be found in 
the later stages of the motion where the vortices shed a I ternati ng I y. 
The accumulation of vorticity reduces the back pressure and gives rise 
to a drag coefficient of about 40 % larger than its steady state value. 
In other words, the initial rise in the drag coefficient is a conse- 
quence of the symmetry of the vortices and the accumulation of 



of vorticity during this symmetric growth. Once the vortices become 
asymmetrical under the influence of ever present small disturbances, 
one vortex moves away from the cylinder causjng a rapid drop In the 
drag force and a gradual increase in the lift force. Evidently the 
rate of rise and the peak value of the drag coefficient are not 
materially affected by the disturbances since the vortices continue 
to retain their symmetry for Ut/c values smaller than about 5. It is 
possible to show that the vortices must acquire certain strength and 
position relative to the cylinder in order to reach a state at which 
they are most susceptible to disturbances. Once the vortices begin to 
shed, the drag coefficient reaches a value of abdut l.l, a value which 
is commonly accepted at the Reynolds numbers encountered herein. 

The peak value of the drag coefficient, its time of occurrence, 
and subsequent decrease to steady state values are also related to the 
character of the fluid motion (laminar or turbulent), and hence on the 
dissipation of circulation. It is a well known fact that the fluid 
motion in the near wake is turbulent for Reynolds numbers larger than 
about 200.- The rate of decrease of circulation in each vortex depends 
on the intensity of turbulence and hence on the Reynolds number. 
Numerical analysis through the use of discrete vortices C3] has shown 
that the larger the dissipation, the smaller is the peak value of 
and the more gradual is the decrease of to its steady state value. 
Evidently, a vortex pair of smaller circulation gives rise to a smaller 
drag coefficient since is proportional to the square of F/Uc. The 
experimental results reported herein are in conformity with the fore- 
going conjectures. At higher Reynolds numbers (Re s 33,000) the peak 
value of Cp is about 1.5 and occurs at a Ut/c value closer to 5, in 
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constrast to a C^-value of 1.55 at Ut/c s; 4 for Reynolds numbers of 
about 16,000. The dissipation of vorticity evidently requires that the 
symmetric vortices be fed a relatively longer time period CA(Ut/c) = 1 3 
in order to aTlow the vortices to reach their maximum strength. As 
noted earlier the vortices gradually become asymmetrical and a turbulent 
vortex, being Inherently more stable to disturbances relative to the one 
at lower Reynolds numbers, lingers somewhat longer in the vicinity of 
the cylinder. These ideas have previously been substantiated through 
measurements and flow visualization by Sarpkaya C33. The extended stay 
of the vortices in the vicinity of the cylinder does not allow the drag 
force to drop rapidly to its steady state value. Suffice it to note 
that the evolution of the drag force with time depends on the Reynolds 
number. The larger the Reynolds number, the larger is the dissipation 
of vorticity, and the smaller the peak value of C^. It has been noted 
previously that there is a plateau in for Ut/c values between 7 and 
13 for larger Reynolds numbers. It is now clear that that plateau is 
a consequence of the dissipation of circulation in the vortices. It is 
also clear that one must further substantiate these ideas by carrying 
out experiments at even higher Reynolds numbers. 

B. LIFT COEFFICIENT 

The lift coefficient versus Ut/c plots are presented in Figs. 9 
through 12 for various values of the Reynolds number. The Reynolds 
numbers in these plots do not exactly correspond to those in the drag 
coefficient plots because of the fact that the lift and drag force 
runs were carried out independently. It was not possible to record 
simultaneously the elevation and the lift and drag forces because of 
the lack of a high-speed three-channel recorder. 
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Figure 9. C vs Ut/c for Reynolds numbers 20,600, 22,100, and 21,700 (two runs). 
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vs Ut/c for Reynolds number 32,600 (two runs). 
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vs Ut/c for Reynolds number 34,100 (three runs). 
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Figure 12. C, vs Ut/c for Reynolds numbers 34,800 and 34,900. 
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A quick examination of the lift plots show that the lift force re- 
mains practically zero until the vortices become asymmetrical at Ut/c 
values of about 6, the amplitude of the lift force varies significantly 
from one run to another even for identical Reynolds numbers, and the 
first few peaks of the alternating lift force do not occur at the same 
Ut/c values. In general, however, the lift coefficient reaches a maxi- 
mum value of about 0.4. This value is in conformity with those obtained 
previously in steady flows at the corresponding Reynolds numbers Cl43. 

There are various reasons for the observed behavior of the lift 
force. Firstly, even in steady flows the lift coefficient shows large 
variations at a given Reynolds number for all Reynolds numbers. This 
is primarily because of the inherent instability of the vortices, the 
surface condition and I ength-to-d i ameter ratio of the cylinder, and 
the lack of spanwise coherence of the vortex cores. Secondly, in im- 
pulsive flow the development of lift is strongly dependent on the 
initial disturbances which give rise to asymmetry and subsequent vortex 
shedding. Consequently, it is not too surprising to find that in some 
runs the vortices continue to grow symmetrically and thus result in 
very small lift forces. In some runs asymmetry develops rapidly and 
the lift force not only starts at smaller values of Ut/c but also 
reaches larger amplitudes. Since there is no way to control the magni- 
tude of the d i sturbances, other than assuring that the fluid in the 
tunnel is as quiescent as possible prior to each run, one must antici- 
pate relatively large variations both in the magnitude and the time 
of occurrence of the peak values of lift. Suffice it to note that the 
lift force is a fairly random quantity and that randomness is certainly 
reflected in the data presented herein. The data further helps to 



explain the random nature of the side forces measured on missile 
models mounted in wind tunnels. It appears that attempts to obtain 
a unique lift-force variation will not be fruitful in view of the 
stated randomness of the lift forces. 
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V. CONCLUSIONS 



The experimental investigation of the impulsively-started flow 
about smooth circular cylinders in a vertical water tunnel warranted 
the following conclusions; 

1. The drag coefficient reaches a peak value of 1.5 1o I. 55, 
depending on the Reynolds number, at Ut/c values from about 4 to 5; 

2. the drag coefficient eventually drops to its steady state 
value for Ut/c larger than about 7; 

3. the lift coefficient begins to develop at Ut/c values larger 
than about 6 and alternates about a zero mean value with an amplitude 
of about 0.4. The magnitude and the time of occurrence of the peak 
values are strongly dependent on the magnitude and distribution of 
the initial disturbances; 

4. the lift and drag coefficients both appear to depend on the 
Reynolds number. Only additional experiments at higher Reynolds 
numbers could determine the extent and form of the dependence on the 
Reynolds number. 
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